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We try a global fit of the experimental branching ratios and CP-asymmetries of the charmless 
B —> PV decays according to QCD factorisation. We find it impossible to reach a satisfactory 
agreement, the confidence level (CL) of the best fit is smaller than .1 %. The main reason 
for this failure is the difficulty to accomodate several large experimental branching ratios of 
the strange channels. Furthermore, experiment was not able to exclude a large direct CP 
asymmetry in — ^ P^i"", which is predicted very small by QCD factorisation. Trying a fit 
with QCD factorisation complemented by a charming-penguin inspired model we reach a best 
fit which is not excluded by experiment (CL of about 8 %) but is not fully convincing. These 
negative results must be tempered by the remark that some of the experimental data used are 
recent and might still evolve significantly. 
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I. INTRODUCTION 



It is an important theoretical challenge to master the non-leptonic decay amplitudes and particularly B non- 
leptonic decay. It is not only important per se, in view of the many experimental branching ratios which have 
been measured recently with increasing accuracy by BaBar [1-10], Belle [11-15] and CLEO [16 21], but it is also 
necessary in order to get control over the measurement of CP violating parameters and particularly the so-called 
angle a of the unitarity triangle. It is well known that extracting a from measured indirect CP asymmetries 
needs a sufficient control of the relative size of the so-called tree (T) and penguins (P) amplitudes. 

However the theory of non-leptonic weak decays is a difficult issue. Lattice QCD gives predictions for semi- 
leptonic or purely leptonic decays but not directly for non-leptonic ones. Since long, one has used what is now 
called "naive factorization" which replaces the matrix element of a four-fermion operator in a heavy-quark decay 
by the product of the matrix elements of two currents, one semi-leptonic matrix element and one purely leptonic. 
For long it was noticed that naive factorization did provide reasonable results although it was impossible to 
derive it rigorously from QCD except in the Nc ^ oo limit. It was also well-known that the matrix elements 
computed via naive factorization have a wrong anomalous dimension. 

Recently an important theoretical progress has been performed [22,23] which is commonly called "QCD 
factorisation". It is based on the fact that the b quark is heavy compared to the intrinsic scale of strong 
interactions. This allows to deduce that non-leptonic decay amplitudes in the heavy-quark limit have a simple 
structure. It implies that corrections termed "non-factorizable" , which were thought to be intractable, can be 
calculated rigorously. The anomalous dimension of the matrix elements is now correct to the order at which 
the calculation is performed. Unluckily the subleading 0{A/mb) contributions cannot in general be computed 
rigorously because of infrared singularities, and some of these which arc chirally enhanced are not small, of 
order O (rn^ / rnt,{rnu + rrid)), which shows that the inverse mt, power is compensated by m^/ (m„ -|- nid)- In the 
seminal papers [22,23], these contributions are simply bounded according to a qualitative argument which could 
as well justify a significanlty larger bound with the risk of seeing these unpredictable terms become dominant. 
It is then of utmost importance to check experimentally QCD factorisation. 

Since a few years it has been applied to i? ^ PP (two charmless pseudoscalar mesons) decays. The general 
feature is that the decay to non-strange final states is predicted slightly larger than experiment while the decay 
to strange final states is significantly Tmd(T(^stimated. In [23] it is claimed that this can be cured by a value of 
the unitarity-trianglc angle 7 larger than generally expected, larger maybe than 90 degrees. Taking also into 
account various uncertainties the authors conclude positively as for the agreement of QCD factorisation with 
the data. In [24,25] it was objected that the large branching ratios for strange channels argued in favor of the 
presence of a specific non perturbative contribution called "charming penguins" [25-30] . We will return to this 
approach later. 

The B PV (charmless pseudoscalar -I- vector mesons) channels are more numerous and allow a more 
extensive check. In ref. [31] it was shown that naive factorisation implied a rather small ]P]/]T] ratio, for 
fiO ^ p^TT^ decay channel, to be compared to the larger one for the B tt+tt". This prediction is still valid 
in QCD factorisation whercjhe \P\/\T\ ratio is of about 3 % (8 %) for the W p+n- P"7r+) chamiel 

against about 20 % for the B° — » tt+tt" one. If this prediction was reliable it would put the B^ p+7r~ channel 
in a good position to measure the CKM angle a via indirect CP violation. This remark triggered the present 
work: we wanted to check QCD factorisation in the B — > PV sector to estimate the chances for a relatively 
easy determination of the angle a. 

The non-charmed B PV amplitudes have been computed in naive factorisation [32], in some extension 
of naive factorisation including strong phases [33] , in QCD factorisation [34-36] and some of them in the so- 
called perturbative QCD [38,39]. In [41], a global fit to B ^ PP,PV,VV was investigated using QCDF in 
the heavy quark limit and it has been found a plausible set of soft QCD parameters that apart from three 
pseudoscalar vector channels, fit well the experimental branching ratios. Recently [36] it was claimed from a 
global fit to _B ^ PP, PV that the predictions of QCD factorization are in good agreement with experiment 
when one excludes some channels from the global fit. When this paper appeared we had been for some time 
considering this question and our feeling was significantly less optimistic. This difference shows that the matter 
is far from trivial mainly because experimental uncertainties can still be open to some discussion. We woiild like 
in this paper to understand better the origin of the difference between our unpublished conclusion and the one 
presented in [36] and try to settle the present status of the comparison of QCD factorisation with experiment. 

One general remark about QCD factorisation is that it yields predictions which do not differ so much from 
naive factorisation ones. This is expected since QCD factorisation makes a perturbative expansion the zeroth 
order of which being naive factorisation. As a consequence, QCD factorisation predicts very small direct CP 
violation in the non-strange channels. Naive factorisation predicts vanishing direct CP violation. Indeed, direct 



CP violation needs the occurence of two distinct strong contributions with a strong phase between them. It 
vanishes when the subdominant strong contribution vanishes and also when the relative strong phase does as is 
the case in naive factorisation. In the case of non-strange decays, the penguin (P) and tree (T) contributions 
being at the same order in Cabibbo angle, the penguin is strongly suppressed because the Wilson coefficients are 
suppressed by at least one power of the strong coupling constant , and the strong phase in QCD factorisation 
is generated by a 0{as) corrections. Having both P/T and the strong phase small, the direct CP asymetries are 
doubly suppressed. Therefore a sizable experimental direct CP asymetry in p+7r~ which is not excluded 

by experiment [9] would be at odds with QCD factorisation. We will discuss this later on. Notice that this 
argument is independent of the value of the unitarity angle 7, contrarily to arguments based on the value of 
some branching ratios which depend on 7 [23]. 

The Perturbative QCD (PQCD) predicts larger direct CP asymmetries than QCDF due to the fact that 
penguin contributions to anihilation diagrams, claimed to be calculable in PQCD, contribute to a larger amount 
to the amplitude and have a large strong phase. In fact, in PQCD, this penguin anihilation diagram is claimed to 
be of the same order, 0{ag), than the dominant naive factorisation diagram while in QCDF it is also ©(a^) but 
smaller than the dominant naive factorisation which is 0(1). Hence, in PQCD, this large penguin contribution 
with a large strong phase yields a large CP asymmetry [40,42,43]. 

If QCD factorisation is concluded to be unable to describe the present data satisfactorily, while there is 
to our knowledge no theoretical argument against it, we have to incriminate non-perturbativc contributions 
which are larger than expected. One could simply enlarge the allowed bound for those contributions which are 
formally subleading but might be large. However a simple factor two on these bounds makes these unpredictable 
contributions comparable in size with the predictable ones, if not larger. This spoils the predictivity of the whole 
program. 

A second line is to make some model about the non-perturbative contribution. The "charming penguin" 
approach [27,30] starts from noticing the underestimate of strange-channels branching ratios by the factorisation 
approaches. This will be shown to apply to the PV channels as well as to the PP ones. This has triggered us to 
try a charming-penguin inspired approach. It is assumed that some hadronic contribution to the penguin loop 
is non-perturbative. In other words that weak interactions create a charm-anticharm intermediate state which 
turns into non-charmed final states by strong rescattering. In order to make the model as predictive as possible 
we will use not more than two unkown complex number and use flavor symmetry in strong rescattering. 

In section II we will recall the weak-interaction effective Hamiltonian. In section III we will recall QCD 

factorisation. In section IV we will compare QCD factorisation with experimental branching ratios and direct 
CP asymmetries. In section V we will propose a model for non-perturbative contribution and compare it to 
experiment. We will then conclude. 



II. THE EFFECTIVE HAMILTONIAN 



The efi^ective weak Hamiltonian for charmless hadronic B decays consists of a sum of local operators Qi 
multiplied by short-distance coefficients Ci given in table I, and products of elements of the quark mixing 
matrix, Xp — VpbV*g or = VpbV*^. Below we will focus on B — > PV decays; where P and V hold for 
pseudoscalar and vector mesons respectively. Using the unitarity relation —At = A„ -|- Ac, we write 

^«ff = % (^1^1 + ^2 Yl CiQi + Cr^Qr^ + CsgQ8g)+h.c., (1) 

p=u,c ^ 1=3,. ..,10 ^ 

where 2 are the left-handed current-current operators arising from PF-boson exchange, Qs^...^ and (57,...,io 
are QCD and elcctroweak penguin operators, and Qr-y and Qsg are the electromagnetic and chromomagnetic 
dipole operators. They are given by 



Ql = 


ipb)v- 


-a{sp)v-a , 




= ipibj)v-A{SjPi)v-A , 


Qz = 


{sb)v- 


-A- {m)v-A , 


Qi 


= {sibj)v-A {qjqi)v-A , 


Q5 = 


{sb)v- 


-A {m)v+A , 


Qe 


= {Sibj)v-A {QjQi)v+A , 


Q7 = 


{sb)v- 


-A^q^qigqW+A, 


Qs 


= {sibj)v-A'Yq ^eg{qjqi)v+A 


Q9 = 


{sb)v- 


-A ^q ^eg{qq)v-A , 


Qw 


= {sibj)v-A'Yg ^eg{qjqi)v-A 


Qr-i = 


^TOf, sa ^,41 +'y5) pi"' b, 


Qsg = 


= ^mf,sa^.(l + 75)G''''6, 



where {qiq2)v±A = 9i7/i(l± 75)92, hj arc colour indices, Cq are the electric charges of the quarks in units of |e|, 
and a summation over q = u,d, s, c, h is implied. The definition of the dipole operators Qj^ and Qsg corresponds 
to the sign convention iD^ = id^ + gsA^ta for the gauge- covariant derivative. The Wilson coefficients are 
calculated at a high scale n ~ Mw and evolved down to a characteristic scale ^ ^ mb using next-to-leading 
order renormalization-group equations. The essential problem obstructing the calculation of non-leptonic decay 
amplitudes resides in the evaluation of the hadronic matrix elements of the local operators contained in the 
effective Hamiltonian. 



NLO 


Ci 


Co 




Ca 




Cfi 


/* — 7T1 1. /'? 
fj/ iit'O / ^ 


1.137 


—0.295 


0.021 


—0.051 


0.010 


—0.065 


H = mb 


1.081 


-0.190 


0.014 


-0.036 


0.009 


-0.042 


H = 2mb 


1.045 


-0.113 


0.009 


-0.025 


0.007 


-0.027 




Ct/q 


Cg/a 


Cg/a 


Cio/a 






jj, = mb/2 


-0.024 


0.096 


-1.325 


0.331 






At = mt 


-0.011 


0.060 


-1.254 


0.223 






/i — 2mj, 


0.011 


0.039 


-1.195 


0.144 






LO 


Ci 


C2 


Ca 


C4 


Cs 


Ce 


H = mb/2 


1.185 


-0.387 


0.018 


-0.038 


0.010 


-0.053 


ti = mb 


1.117 


-0.268 


0.012 


-0.027 


0.008 


-0.034 


yu = 2mb 


1.074 


-0.181 


0.008 


-0.019 


0.006 


-0.022 




Ct/q 


Cs/a 


Cs/q 


Cio/a 






H = mb/2 


-0.012 


0.045 


-1.358 


0.418 


-0.364 


-0.169 


^l = mb 


-0.001 


0.029 


-1.276 


0.288 


-0.318 


-0.151 


fi = '2iiii, 


0.018 


0.01!) 


-1.212 


0.1!):! 


-0.281 


-0.136 



TABLE I. Wilson coefficients d in the NDR scheme. Input parameters are = 0.225 GeV, mt{mt) = 167 GeV, 
mb{mb) = 4.2 GeV, Mw = 80.4 GeV, a = 1/129, and sm'^9w = 0.23 [23]. 



III. QCD FACTORIZATION IN B ^ PV DECAYS 



When the QCD factorization (QCDF) method is appUed to the decays B—^PV, the hadronic matrix elements 
of the local effective operators can be written as 

( PV\Oi\B) = Ff^PiO)Tl^,* fv^v + A^^"" iO)Tj,^i* fp^P 

+ Tl'*fB<^B*fv^v*fp<^p, (3) 



where $m arc leading-twist light-cone distribution amplitudes, and the ^-products imply an integration over 
the light-cone momentum fractions of the constituent quarks inside the mesons. A graphical representation of 
this result is shown in Figure 1. 

Here F^^^ and Aq^^ denote the form factors for B^P and B^V transitions, respectively. $b(C), $y(a;), 
and <i>p(y) are the light-cone distribution amplitudes (LCDA) of valence quark Pock states for B, vector, and 
pseudoscalar mesons, respectively. T"/'^^ denote the hard-scattering kernels, which are dominated by hard gluon 
exchange when the power suppressed 0{Aqcd/ mb) terms are neglected. So they are calculable order by order in 
perturbation theory. The leading terms of come from the tree level and correspond to the naive factorization 
(NF) approximation. The order of terms of can be depicted by vertex-correction diagrams Fig. 2 (a-d) 
and penguin-correction diagrams Fig. 2 (e-f). T^^ describes the hard interactions between the spectator quark 
and the emitted meson M2 when the gluon virtuality is large. Its lowest order terms are 0{as) and can be 
depicted by hard spectator scattering diagrams Fig. 2 (g-h). One of the most interesting results of the QCDF 
approach is that, in the heavy quark limit, the strong phases arise naturally from the hard-scattering kernels 
at the order of ag. As for the nonperturbative part, they are, as already mentioned, taken into account by the 
form factors and the LCDA of mesons up to corrections which are power suppressed in l/mi,. 




With the above discussions on the effective Hamiltonian of B decays Eq.(l) and the QCDF expressions of 
hadronic matrix elements Eq.(3), the decay amplitudes for B^P V in the heavy quark limit can be written as 

„ 10 

A{B^pv) = E vmB)^^. (4) 

* p=u,c i=l 

The above {PV\Oi\B)^j: are the factorized hadronic matrix elements, which have the same definitions as those in 
the NF approach. The "nonfactorizable" effects are included in the coefficients ai which are process dependent. 
The coefficients are collected in Sec. Ill A, and the explicit expressions for the decay amplitudes of B^P V 
can be found in the appendix A. 

According to the arguments in [22], the contributions of weak annihilation to the decay amplitudes are 

power suppressed, and they do not appear in the QCDF formula Eq.(3). But, as emphasized in [40,42,43], the 
contributions from weak annihilation could give large strong phases with QCD corrections, and hence large CP 
violation could be expected, so their effects cannot simply be neglected. However, in the QCDF method, the 
annihilation topologies (see Fig. 3) violate factorization because of the cndpoint divergence. There is similar 
endpoint divergence when considering the chirally enhanced hard spectator scattering. One possible way is 
to treat the endpoint divergence from different sources as different phenomenological parameters [23]. The 
corresponding price is the introduction of model dependence and extra numerical uncertainties in the decay 
amplitudes. In this work, we will follow the treatment of Ref. [23] and express the weak annihilation topological 
decay amplitudes as 



A'^iB^PV) (X Jb fp fvJ^Xp , 



(5) 



where the parameters bi are collected in Sec. Ill B, and the expressions for the weak annihilation decay ampli- 
tudes of B^P V are listed in the appendix B. 



A. The QCD coefficients Oi 



We express the QCD coefficients Oj (see Eq.(4) ) in two parts, i.e., Ui = Uij + aijj. The first term Uij 
contains the naive factorisation and the vertex corrections which are described by Fig. 2 (a-f), while the second 
part Uiji corresponds to the hard spectator scattering diagrams Fig. 2 (g-h). 

There are two different cases according to the final states. Case I is that the recoiled meson Mi is a vector 
meson, and the emitted meson M2 corresponds to a pseudoscalar meson, and vice versa for case II. For case I, 
we sum up the results for Ui as follows: 
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1 + 
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An 

CpOts 
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+ 



9n Nc ' 



Vm 



Nc 



1 + ^Vm 
An 



Qtt Nc 



02,77 = 
03,77 = 
04,77 
05,77 = 



iV2 

nCpUg 
nCpfXs 

~w 

nCpct. 
nCpoc 



N} 



CiH{BMi,M2), 
CiH{BMi,M2), 
'-C3H{BMi,M2), 
CeH'{BMi,M2), 



0.6,11 = 0, 



nC 



07,77 = -- 



pQ-s 



asji = 0, 

nCpcts 



09,77 



010,77 



7V2 

nCpoi 



C8H'{BMi,M2), 

CwH{BMi,M2), 
'-CsH{BMi,M2), 



(6) 



where Cp = %v~' ^^'^ ~ "^^^ vertex parameters Vm and result from Fig. 2 (a-d); the QCD penguin 
parameters ^ and the electroweak penguin parameters result from Fig. 2 (e-f). 

The vertex corrections are given by: 



Vm 


= 121n^- 




M 


V'm 


= 121nl^- 




9{x) 







2Li2(a;) - In^a; 



2 In a; 



(3 + 2in) Inx — [x \ — x) 



(7) 



where \a2[x) is the dilogarithm function, whereas the constants 18 and 6 are specific to the NDR scheme. 
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FIG. 2. Order corrections to the hard-scattering kernels. The upward quark hues represent the emitted mesons from 
b quark weak decays. These diagrams are commonly called vertex corrections, penguin corrections, and hard spectator 
scattering diagrams for Fig. (a-d), Fig. (e-f), and Fig. (g-h) respectively. 



The penguin contributions are: 

[ 



pp 

M,2 



pP 



+ X] + Ce + ^eqCs - 

q=q' 

^ [4^ mi, 2 - , , 
Ci[-ln-^ + --Gm(sp) 



C3 — 



^ln^ + |-GM(0)-GM(l) 



-CqClQ 



4, TO6 

^In GnKSq) 



-2^8^/ dx 



1-a; 



1 



C-^ — — Cq 

J 2 



^ln^ + ^-GM(0)-GM(l) 



- ^ (g4 + Ge + ^e,G8 + ^e,Gi 

q=q' 



GAf(sq) 



(8) 



and the electroweak penguin parameters P^f^ '■ 



= {c, + N,c,) 



^In— + ^ - Gm(sp)] - (G3 + N,C, 



3 



|n!^ + ^4G..(0)4GM(l) 



^ (tV.Gs + G4 + A^cGs + Ge) -e, [-In^ - GM(sg)] - N,C',^ / '^^ T 



rd mi 2 - 1 I 

In— + ^ - Gm(sp)J - (G3 + iV,G4 



r^,.I^ + |-iG„,o)-ieM(i) 



+ (tVcGs + G4 + A^cGs + G 



3 [4 mb 



In— - GA/(sg) 
3 ^ 



(9) 



where Sq — / , and where in the expressions for ^ and P^f'^ runs over all the active quarks at the 
scale /i = 0{mi,), i.e., q' — u, d, s, c, 6. The functions Gm{s) and Gm(s) are given respectively by: 



Gm(s) = 


[ dxG{s~ 




Jo 






Gm(s) = 


[ dxG{s^ 




Jo 






G{s,x) = 


—4 / duu 




Jo 



-My 



2(12s + 5a;-3a;lns) 4V4s -x{2s + x) 



9x 



3x3/2 



arctan 



4s — X 



(10) 
(11) 

(12) 



The parameters H{BMi, M2) and H'{BMi, M2) in aijj, which originate from hard gluon exchanges between 
the spectator quark and the emitted meson M2, are written as: 



H{BV, P) 



fB.h 



m 



H'{BV,P) = 



Jo Jo 

dx dy- 

{0)J J J ? X y 



^ X y '' 
$b(0 ^p{x) ^y(2/) 



(13) 



For case II (vector meson emitted) except for the parameters of H{BMi, M2) and H'{BMi, M2), the expres- 
sions for a, are similar to those in case I. However we would like to point out that, because {V\{qq)s±p\0) = 0, 
the contributions of the effective operators Oe.s to the hadronic matrix elements vanish, i.e., the terms that 
are related to 05,8 disappear from the decay amplitudes for case II. As to the parameters H{BMi,M2) and 
H'{BMi,M2) in o,,//, they are defined as 



H{BP, V) 
H'{BP, V) 



IbIf 



m|Ff-^(0) 

fsfp 
m|Ff-^(0) 



1 .1 
dx / dy 

^0 
1 rl 



$s(0 ^v{x) \ ^p{y) ^ 2/ip X $^(y)] 



rrib x y 



dS, dx dy 



e X 

$b(0 <^v{x) r$p(y) ^ 2np x $P (y) 



y 



+ ■ 



rub x y 



(14) 



The parameter = m'j,/{mi + 1112) where mi, 2 are the current quark masses of the meson constituents, is 
proportional the the chiral quark condensate. 



B. The annihilation parameters bi 



The parameters of bi in Eq.(5) correspond to weak annihilation contributions. Now we give their expressions, 
which are analogous to those in [23]: 

6i(Mi,M2) = ^CiA\{Mi,M2), 

62(Mi,M2) = ^C2A\{M,,M2), 



bs{M„M2) 
bi{MuM2) 



Cf 
Cf 



^2 { CsAl (Ml , M2) + C5 A^(Mi , M2) + [C5 + N,Ce\ A^ (Mi , M2) } , 



^2 {C4^i(Mi, M2) + ^^(Mi, M2)}, 
bl^{M,,M2) = ^{CgAUMi, M2) + C7A^(Mi, M2) + [Cr + iV,C8]4(Mi, M2)}, 
6r(Mi,M2) = ^{Cio4(Mi,M2) + C84(Mi,M2)}. 



(15) 



Here the current-current annihilation parameters 61,2(^^1, M2) arise from the hadronic matrix elements of the 
effective operators ©1,2, the QCD penguin annihilation parameters 63,4(Mi, M2) from O3-Q, and the electroweak 
penguin annihilation parameters 63'4(Mi, M2) from Ot-iq- The parameters of bi are closely related to the final 
states; they can also be divided into two different cases according to the final states. Case I is that Mi is a 
vector meson, and M2 is a pseudoscalar meson (here Mi and M2 are tagged in Fig. 3). Case II is that Mi 
corresponds to a pseudoscalar meson, and M2 corresponds to a vector meson. For case I, the definitions of 
4'^(Mi,M2) in Eq.(15) are 



AyV,P) = 0, 

Ai{V,P) = TraJ dxf dy<^vix)'S>Uy) 



2/ip 2{l + x) 



Jo Jo 

A\{V,P) = na, [dxf dy^v{x)^p{y) 
Jo Jo 

Ai{V, P) = -TTUs f dxf dy<Pv{x)<^p{y) 
Jo Jo 

Ai{V,P) = nas f'dx f\y^v{xWp{y) 



pyyj -2 
mb x^y 



+ 



y{l - xy) x'^y 



+ 



1 



Jo 



x{l — xy) xy"^ 
2fJ.p 2y 
rub xy{l - xy) ' 



(16) 



For case-II, 



Jo Jo 



A\{P, V) = nas f dx [ dy$p(x)$y(y) 
Jo Jo 

Ai{P, V) = -nas ( dx I dy^p{x)^v{y) 
Jo Jo 

AliP, V) = TTUs [ dx [ dy^Pp{x)<S>v{y) 
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rrib xy^ 
1 1 



y{l - xy) x'^y 
1 1 
x{l — xy) xy"^ 
2^,p 2x 
nib xy{l - xy) 



(17) 



Here our notation and convention are the same as those in [23]. The superscripts i and / on A^'^ correspond 
to the contributions from Fig. 3(a-b) and Fig.3(c-d), respectively. The subscripts k — 1,2,3 on A^i^^ refer to 
the Dirac structures {V - A)®{V - A), {V ~ A)®{V + A) and {-2){S - P)®{S + P), respectively, ^vix) 
denotes the leading- twist LCDAs of a vector meson, and ^p{x) and ^^p{x) denote twist-2 and twist-3 LCDAs 
of a pseudoscalar meson, respectively. 




laj (b) (c) (d) 

FIG. 3. Order as corrections to the weak annihilations of charmless decays B~ 



>PV. 



Note that assuming SU{3) flavor symmetry implies symmetric LCDAs of light mesons (under x-^x), whence 
A\ — —A\. In this approximation the weak annihilation contributions (for case I) can be parametrized as 



Al(F,F)~187ra,(x^-4+y 



M,{V,P)^i:asr^ 



2tt^ -Qixl 



Aily, P)~6nasr^ (2X1 " 



2Xa 



(18) 



where Xa = f^dx/x parametrizes the divergent endpoint integrals and = 2/ip/mb is the so-called chirally 
enhanced factor. We can get similar forms to Eq.(18) for case II, but with A{{P,V) = ~A{{V,P). In our 
calculation, we will treat Xa as a phenomenological parameter, and take the same value for all annihilation 
terms, although this approximation is crude and there is no known physical argument justifying this assumption. 
We shall see below that Xa gives large uncertainties in the theoretical prediction. 



IV. QCD FACTORISATION VERSUS EXPERIMENT 



In order to propose a test of QCD factorization with respect to experiment, a compilation of various charmless 
branching fractions and direct CP asymmetries was performed and is given in tables II, III and IV. This 
compilation includes the latest results from BaBar, Belle and CLEO. The measurements were combined into a 
single central value and error, that may be compared with the theoretical prediction. First, the total error from 
each experiment was computed by summing quadratically the statistic and systematic error: this approach is 
valid in the limit that the systematic error is not so large with respect to the statistic error. Secondly, when the 
experiment provides an asymmetric error i ^ conservative symmetric error was assumed in the calculation 
by using ±Max(cri, CT2)- In case of a disagreement between several experiments for a given measurement, the 



total error was increased by a "scale factor" computed from a combining the various experiments, using the 
standard procedure given by the PDG [44]. 

In order to compare the theoretical predictions {y} with the experimental measurements {xia^}, the following 
was defined: 

In the case when a correlation matrix between several measurements is given by the experiment, as in the case 
of the p'^-K~ / p'^K~ measurements, the was corrected to account for it. The above was then minimized 
using MINUIT [46], letting free all theoretical parameters in their allowed range. The quality of the minimum 
yielded by MINUIT was assessed by replacing it with an ad hoc minimizer scanning the entire parameter space. 
The theoretical predictions, with the theoretical parameters yielding the best fits, are compared to experiment 
in table VI for two scenarios to be explained below. The asymmetries of the p^7r^ channels can be expressed [9] 
in terms of the quantities reported in table IV. The comparison between their theoretical predictions and 
experiment is reported in table VII. 

Scenario 1 refers to a fit according to QCD factorisation, varying all theoretical parameters in the range 
presented in table V. Even the unitarity triangle angle 7 is varied freely and ends up not far from 90 degrees. 
We have taken Xa = Xh in the range proposed in ref. [23] : 

1 

Xa,h = / ^ = In ^(1 + pA,H e'^-'-). (19) 



These parameters label our ignorance of the non perturbatively calculable subdominant contribution to the 
annihilation and hard scattering, defined in Eqs. (16, 17) and Eqs. (13, 14) respectively. They do not need to 
have the same value for all PV channels but we have nevertheless assumed one common value since a fit would 
become impossible with too many unknown parameters. 

Scenario 2 in table VI refers to a fit adding a charming penguin inspired long distance contribution which will 
be presented and discussed in section V. In this fit 7 is constrained within the range [34°, 82°]. 

The values of the theoretical parameters found for the two best fits is given in table V: many parameters 
are found to be at the edge of their allowed range^^. In order to estimate the quality of the agreement between 
measurements and predictions, the standard Monte Carlo based "goodness of fit" test was performed: 

• the best-fit values of the theoretical parameters were used to make predictions for the branching ratios 
and CP asymmetries, 

• the total experimental error from each measurement was used to generate new experimental values dis- 
tributed around the predictions with a Gaussian probability, 

• the full fit previously performed on real measurements is now run on this simulated data, and the of 
this fit is saved in a histogram H. 

It is then possible to compare the Xdata obtained from the measurement with the x^ one would obtain if the 
predictions were true. Additionally, one may compute the confidence level of the tested model by using: 

The results of the "goodness of fit" tests are given in FIG. IV. From these tests, one may quote an upper 
limit for the confidence level in scenario 1, CL < 0.1%, and in the case of scenario 2, CL < 7.7%. 



^^Tablc 5 shows that the fit value of pA appears at the edge of the input range, pA = 1- However enlarging the range 
of Pa, such as \pa\ < 10, brings a large annihilation contributions {pA, 4>a) = (2.3, —41°) for scenario 1 and (4.4, —108°) 
for scenario 2. With so large values of \pa\ the unpredictable contributions would dominate the total result making the 
whole exercize void of signification. 




(a) Scenario 1 (b) Scenario 2 

FIG. 4. Goodness of fit test of the two proposed theoretical models: the arrow points at the value Xdata found from 
the measurements, and the histogram shows the values allowed for x'^ in the case that the models predictions are correct. 
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TABLE II. Experimentally known data of CP-averaged branching ratios for the charmless B PV decay modes, 
used as input for the global fit. The channels containing the r)' meson have been excluded. 
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TABLE III. Experimental measured data of direct CP asymmetries for the charmless B — > PV decay modes, used as 
input for the global fit. 
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TABLE IV. Experimental results and correlation matrix for the various asymmetries measured in the channels 
p^-K^ / p^K^ . The notations are explained in [9]. 



In tables II (III) we give the experimental CP-averaged branching ratios (direct CP asymmetries) which we 
have used in our fits. We have also used the quantities reported in table IV which are related to the branching 
ratios and CP asymmetries of the B — > p^-K^ channels. 
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TABLE V. Various theoretical inputs used in our global analysis of i? — > PV decays in QCDF. The parameter ranges 
have been taken from literature [23], [34], [35] and [37]. The two last columns give the best fits of both scenarios. 
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TABLE VI. Best fit values using the global analysis of B ^ PV decays in QCDF with free 7 (scenario 1) and 
QCDF±Charmirig Penguins (scenario 2) with constrained 7. The CP-averaged branching ratios are in unit of 10^®. 





Experiment 


Scenario 1 


Scenario 2 


,p+ IT 


-0.82 ±0.31 ±0.16 


-0.04 


-0.23 


AP~ 


-0.11 ±0.16 ±0.09 


-0.0002 


0.04 



TABLE VIL Values of the CP asymmetries for B ^ wp decays in QCDF (scenario 1) and QCDF±Charming Penguins 
(scenario 2). The notations are explained in [9]. 



For the sake of definiteness let us remind that the branching ratios for any charmless B decays, B — > PV 
channel, in the rest frame of the B-meson, is given by 

BTZ{B ^ PV) = '^-^\A{B ^ PV)+A"-{B ^ PV)+A^^{B ^ PV)\'^, (20) 

where tb represents the i3- meson life time (charged or uncharged according to the case). The amplitudes 
A, A"' and A^^ are defined in appendix A, B and in eqs. (24) and (25) respectively. In the case of pure QCD 
factorisation (scenario 1) we take of course A^^ = 0. The kinematical factor |p| is written as: 

I I ^ v/[m^ - (mp + mvy][mjj - {mp - 



A. Comparison with Du et al 



Our negative conclusion about the QCD factorisation fit of the B — > PV channels is at odds with the 
conclusion of the authors of ref. [36], who have performed a successful fit of both B PP and B PV 
channels using the same theoretical starting point. These authors have excluded from their fits the channels 
containing a K* in the final state, arguing that these channels seemed questionable to them. We have thus 
made a fit without the channels containing the K*, and indeed we find as the authors of ref. [36] that the global 
agreement between QCD factorisation and experiment was satisfactory. Notice that this fit was done without 
discarding the channels — > u!tt'^{K~^) as done by Du et al. 

Notice also that the parameters Cpjr and the A(^p have been kept in this fit. The disagreement between 
QCDF and experiment for these quantities was not enough to spoil the satisfactory agreement of the global fit 
because the experimental errors arc still large on these quantities. 

The conclusion of this subsection is that the difference between the "optimistic" conclusion about QCDF of 
Du et al and our rather pessimistic one comes from their choice of discarding the channels containing the K*'s. 
In other words the conclusion about the status of QCDF in the B PV channels depends on the confidence 
we give to the published results on these channels. 



V. A SIMPLE MODEL FOR LONG DISTANCE INTERACTIONS 



As seen in table VI the failure of our overall fit with QCDF can be traced to two main facts. First the 

strange branching ratios are underestimated by QCDF. Second the direct CP asymmetries in the non-strange 
channels might also be underestimated. A priori this could be cured if some non-perturbative mechanism was 
contributing to |P|. Indeed, first, in the strange channels, |P| is Cabibbo enhanced and such a non-perturbative 
contribution could increase the branching ratios, and second, increasing |P|/|T| in the non-strange channels 
with non-small strong phases could increase significantly the direct CP asymmetries as already discussed. We 
have therefore tried a charming-penguin inspired model. We wanted nevertheless to avoid to add too many new 
parameters which would make the fit void of signification. We have therefore tried a model for long distance 
penguin contributions which depends only on two fitted complex numbers. 

Let us start by describing our charming-penguin inspired model for strange final states. In the "charming 
penguin" picture the weak decay of a B {B~) meson through the action of the operator (see notations 

in Eqs. (1) and (2)) creates an hadronic system containing the quarks s,d{u),c,c, for example D^*^ + D^*'> 
systems. This system goes to long distances, the c, c eventually annihilate, a pair of light quarks are created by 
non-perturbative strong interaction and one is left with two light meson. Let us here restrict ourselves to the 
case of a PV pair of mesons, i.e. one of the final mesons is a light pseudoscalar (tt, K, rf) and the other a light 
vector meson (p, tj, 0, K*). We leave aside from now on the r]' which is presumably quite special. 

We will picture now this hadronic system as a coherent state which decays into the two final mesons with 
total strangeness -1. This state has a total angular momentum J = 0. Its flavor sd is that of a member of an 
octet of flavor- S'f/(3). We will assume flavor-5C/(3) symmetry in the decay amplitude of this hadronic state. 
This still leaves four SU (3)-invariant amplitudes since both P and V can have an octet and a singlet component 
and that there exist two octets in the decomposition of 8 x 8. We make a further simplifying assumption based 
on the OZI rule. Let us give an example: we assume that V = {sq) where q is any of the light quarks u, d, s, 
and that P = (qd). Then we compute the contractions between 



< {sq){qd)\s{uu + dd + ss)d >= 1 



(22) 



The meaning of this rule is simple. We add to the sd quarks in our hadronic state an SU (3) singlet uu+dd+ss 
and compute an "overlap" making contractions so that the quarks in the singlet go into two different mesons. 
This latter constraint is the OZI rule. This is why the overlap in Eq. (22) is 1 even \i q = d since it is forbidden 

to have both d quarks from the singlet in the same final meson. As an example, the decay B ^ K gives the 
following overlap coeflacient: 

, ^, (uu — dd) , , - , , - 1 , 

<{sd)- -j^-^\s{uu + dd + ss)d >= p (23) 

V 2 v2 

For the r] meson we will use the decomposition in [32] . The overlap coefficients thus computed play the role 

of SU{3) Clebsch-Gordan (CG) coefficients computed in a simple way. These coefficients are assumed to be 
multiplied by an universal complex amplitude to be fitted from experiment. Up to now we have assumed that 
the active quark (here s) ended up in the vector meson. We need another universal amplitude for the case where 
the active quark ends up in the pseudoscalar meson. 

We are thus left with two theoretically independent and unknown amplitudes, one with V — (sq), P — (qd), 
one with P — (sq), V = (qd). We shall write them respectively as (-4^) when the active quark ends up in 
the Pseudoscalar (Vector) meson. 

Concerning the B decay into a pseudoscalar + vector meson of vanishing total strangeness, we apply the 
same recipe with the same amplitudes and X^, replacing the s quark by a d quark and, of course, the 
corresponding replacement of the CKM factor VcbV*^ by VcbV*^. 

a. To summarize : the long distance term is given by two universal complex amplitudes multiplied by a 
CG coefiicient computed simply by the overlap factor in (23), see table VIII. 
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TABLE VIII. F\avoT-SU (3) Clebsch-Gordan coefficient for long distance penguin-like contributions. Notice that the 
channel B~ — » $ 7r~ vanishes due to OZI rule. 



In practice, to the amplitudes described in the appendices we add the long distance amplitudes, given by: 



A^^{B ^ PV) = ^ml\'^{ClP A"" + Cl^ A^) (24) 
for the non strange channels and 



^LD(B ^ PV) = ^mlXciClP A^ + Cl^ A^) (25) 
v2 

for the strange channels. In equations (24) and (25), A^ and A^ arc two complex numbers which arc fitted in 
the global fit of scenario 2 and Cl^ and Cl^ are the &a.voi-SU (3) Clebsch-Gordan coefficients which are given 
in table VIII. For both channels containing the r] we have used the formulae 

^ cos^ _ sin^ C/^ = _2£^_!^ (26) 

with 00 = -9.1° and = -22.2°. 

The fit with long distance penguin contributions is presented in tabic VI under the label "Scenario 2" . The 
agreement with experiment is improved, it should be so, but not in such a fully convincing manner. The goodness 
of the fit is about 8 % which implies that this model is not excluded by experiment. However a look at table V 
shows that several fitted parameters are still stuck at the end of the allowed range of variation. In particular 
Pa = ^ means that the uncalculable subleading contribution to QCDF is again stretched to its extreme. 

Finally the fitted complex numbers which fix the size of the long distance penguin contribution (last four lines 
in table V) are small. To make this statement quantitative, assuming the long distance amplitude were alone, 
the values for A^ and A^ in table VI correspond to branching ratios which reach at their maximum 6 x 10~^ 
but are more generally in the vicinity of 2 x 10~^. In part, this is due to the fact that, if some strange channels 
want a large non-perturbative contribution to increase their branching ratios, some other strange channels and 
particularly the B K<j) channels which are in good agreement with QCDF cannot accept the addition of a too 
large non-perturbative penguin contribution. This last point should be stressed: if the strange channels show 
a general tendency to be underestimated by QCDF, there is the striking exception of the sss channels which 
agree very well with QCDF and make the case for charming penguins rather difficult. 



VI. CONCLUSION 



We have made a global fit according to QCD factorisation of published experimental data concerning charmless 
B PV decays including CP asymmetries. We have only excluded from the fit the channels containing the 
rj' meson. Our conclusion is that it is impossible to reach a good fit. As can be seen in the scenario 1 of table 
VI, the reasons of this failure is that the branching ratios for the strange channels are predicted significantly 
smaller than experiment except for the B (f>K channels, and in table VII it can be seen that the direct CP 
asymmetry of _B ^ p+Tr^ is predicted very small while experiment gives it very large but only two sigmas from 
zero. Not only is the "goodness of the fit" smaller than .1 %, but the fitted parameters show a tendency to 
evade the allowed domain of QCD factorisation. One might wonder if wc were not too strict in imposing the 
same scale /i in all terms since the value of /i, representing the effect of unknown higher order corrections, could 
be different in different classes of channels We have performed several tests relaxing this unicity of fi and 
concluded that it affected very little the outcome of our fit. 

For the sake of comparison with the authors of ref. [36] we have tried a fit without the channels containing a 
K* . The result improves significantly. The only lesson we can receive from this is that one must look carefully 

at the evolution of the experimental results, many of them being recent, before drawing a final conclusion. 

Both the small predicted branching ratios of the strange channels and the small predicted direct CP asym- 
metries in the non strange channels could be blamed on too small P amplitudes with too small "strong phases" 
relatively to the T amplitudes. We have therefore tried the addition of two "charming penguin" inspired long 
distance complex amplitudes combined, in order to make the model predictive enough, with exact flavor-^J/ (3) 
and OZI rule. This fit is better than the pure QCDF one: with a goodness of the fit of about 8 % the model 
is not excluded by experiment. But the parameters show again a tendency to reach the limits of the allowed 



We thank Gerhard Buchalla for raising this question. 



domain and the best fit gives rather small value to the long distance contribution. The latter fact is presumably 
due to the B (j)K which are w(^ll predicted by QCDF and thus deliver a message^ which contradicts the other 
strange channels. This seems to be the reason of the moderate success of our "charming penguin" inspired 
model. 

Altogether, the present situation is unpleasant. QCDF seems to be unable to comply to experiment. QCDF 
implemented by an ad-hoc long distance model is not fully convincing. No clear hint for the origin of this 
problem is provided by the total set of experimental data. PQCD, also called kr factorisation, would predict 
larger direct CP asymmetries, but we do not know if their sign would fit experiment neither if an overall 

agreement of the branching ratios with data can be achieved. 

Maybe however, the coming experimental data will move enough to resolve, at least partly, this discrepancy. 
We would like to insist on the crucial importance of direct CP asymmetries in non-strange channels. If they 
confirm the tendency to be large, this would make the case for QCDF really difficult. 

Finally we do not know yet the answer to our initial question: arc wc in a good position to study the unitarity- 
triangle angle a from indirect CP asymmetries thanks to small penguins. If experimental data evolve so as to 
provide a better support to QCDF, one could become bold enough to use it in estimating a and this would 
reduce the errors. Else, only model-independant bounds [45] could be used but they are not very constraining 
in part because of discrete ambiguities. 
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APPENDIX 

APPENDIX A: THE DECAY AMPLITUDES FOR B^PV 

Following ref. [32], we give the decay amplitudes for the following B—^PV decay processes: 
(1) b ^ d processes: 



^ p-7r+) = ^m|/,Ff ^-(m^){A>i + (A^ + X'J[a4 + aio]}- (Al) 

^ p+TT-) = ^m|^<^^(m2){A;ai + (A^ + K)[a4 + aio - r^ae + as)]}. (A2) 

A(B° 7rV°) = [f^A^-'iml) |a; a2 - + K) [«4 - ^aio - r^iae - ^ag) + - ag)] | 

+f^Ff-^ml) |a>2 - (A; + K)[a, - ^aio - ^(ar + ag)]}) • (A3) 
A{B- ^ TT-pO) = [uA^-'imD {A^i + (A^ + X'Jia^ + - r^iae + as)]} 

+ fpFf-^ml) |A;a2 + (A^ + A'J[-a4 + + ^(07 + a9)]|) . (A4) 
A{B- ^ p-n°) = ^m|(^<^^K) |a; a2 + (A^ + A'J -a, + ^a,o - r^i-ae + ^as) + ^(ag - ar) | 

+/pFf (ml){X'^a, + i\'^ + \'J[a, + aio] }) • (A5) 
A{B- ^ TT-io) = ^ml(^UA^^'^{ml) {A>i + {K + K)[ai + a^o - r^ia^ + as)]} 



+ f.F^-''{ml)^\'^a2 + {X'^ + K) 



04 -I- 2(03 -I- as) + ^(a7 -I- 09 - aio) 



(A6) 



(2) 6 — > s processes: 



A{B° ^ K*-w+) = ^mlfK^Fi'^^ml,)[\uai + (A„ + Ae)[a4 + aio]}. 



V2 



A{B° K-p+) = ^m|/if A^^''(m|f ){a„oi + (A„ + Ac)[a4 + aw - (ag + ag)]}. 



A{B° K°p°) = ^rn% |/if^o^^''(TOjs:o)(-A„ - Ac) 04 - ^aio - r^iae - ^ag) 



Ka2 + (A„ + Ac) X i^{aQ + a?) 



A{B- ^ i^*-7r°) 



Gf 



-nil 



UA^^^ (mi) { A„a2 + (A„ + Ac) x -(ag - ar) 



} 



+ fK'F^^^{m\,){\uai + (A„ + Ac)(a4 + aw)}] ■ 



Gf 2 



fKAQ^P{m%) {A„ai + (A„ + Ac)[a4 + oio - r^{ae + as)]} 



+ fpFF^^iml) |A„a2 + (A„ + Ac) x ^(ag + 07) | 



A{B" ^ ITij) = ( /K^^^'^(m^o)(A„ + Ac 



+ /ci^f {mi) <! A„a2 + (A„ + Ac) 
Gf 2 



1 1 N 



2(03 + 05) + -(Og +07) 



-m, 



fKAQ^'^{m%) {A„ai + (A„ + Ac)(a4 + aw - r'^iae + as))} 



{mi) I A„a2 + (A„ + Ac) |^2(a3 + 05) + ^{aa + 07) 
A{B- ^ i^*-77W) = ^ml[fK'F^-^"\ml) {A„ai + (A„ + Ac)(a4 + aio)} + ^0 ) 

2(a3 - 05) + ^(09 - ay) + ^y^'^ae - ;^a8) 



A„a2 + Ka2-^ + {K + Ac) 

/„(') 



(A7) 
(A8) 

(A9) 

(AIO) 
(All) 

(A12) 

(A13) 
(A14) 



+(03 - 05 + ag - 07) 4r" + ('^a - 05 - ^(ag - 07) + 04 - ^aw - 'r'^'\aQ - ^as)") 4r^ 

ry(') N /J 7y(') 



^(# ^ r%W) = (/K.Ff (m^)(A„ + Ac) ^ ' 



04 — 2*^10 



+/;(o^r''*KM) 



(A15) 



A„a2 + Aca2-^ + (A„ + Ac) 
/,,(') 

+(a3 - 05 + ag - 07) + ( as - as 
/,,(') 



2(a3 - as) + ^(ag - 07) + rf'^{a& - \a^)^ 



{az - as - ^(ag - a-j) + 04 - ^aio - r-^J^'^ae - \°'^) \ -pT 



(3) Pure penguin processes: 



A{B- 7r-K*°) = %mlfK'F^^%m],,){\u + \c) 



V2 



0-4 — -^aiQ 



^0. G 



A{B- - K-K*°) = ^m|/K.Ff-^(m^.)(A; + K) 
A{B- K*-K^) = %mlfKA^-^' {mlo){K + K) 



V2 

A{B- ^ n-cf>) = -^mlUF^-^ml){\'^ + A^,; 
A{B- ^K-(fy= A(^ ifcf)) 



1 

04 — 2^10 



013 + as - 2(^7 + 09) 



as + a4 + as - -{a^ + ag + aio) 



APPENDIX B: THE ANNIHILATION AMPLITUDES FOR B^PV 



We give in this section the following annihilation amplitudes for B^PV already given in ref. 

different notations: 



(1) b ^ d processes: 



-^6r(7r-,p+) + 6r(p+,7r-)-i6r(7r-,p+) 



A%B ^TT+p-) = -^/BA/p|A;6i(7r+,p-) + (A; + A;)[63(p-,7r+)+64(7r+,p-) + 64(p 
- \ bT (P- , ) + br (7r+ , p- ) - i 6r (p- , 7r+ ) 



Gf 



fBUfp\K 6i(p",^")+6i(^",p") 



+(A; + a;) [63(p°, 7r°) + 63(7r°, p°) + 264(77°, p°) + 264(p°, 7r°) 

+1 ( - 6r(p°, - bri^", p°) + 6r (tt", p°) + 6r (p°, tt^))] |. 

^''(B-^7r-p°) = ^/B/./p|A;[62(7r-,p°) - 62(p°,7r-) 



+(A; + K) [63(7r-, pO) - 63(p°, 7r-) + 6^ (7r-, p°) - 6^ (p°, 7r- 



^"(B-^ttV) = ^fBfM a; 62(p-,7r°) - 62(77°, p-) 



+ (AU A^) 63(p-,7rO) - 63(^0, p-) + br{p-,n") - br{7r",p-) 



A''{B-^7T-w) = ^feUUi K 62(77-, w) +62(0;, 7r- 



+ (A; + A;) 63(^-,u;)+63(u;,7r-)+6r(7r-,u;) + 6r(u;,^-) \. 



(2) 6 — > s processes: 



A''(B'^K-p+) = %fsfKfp\ (A„ + Ac) h^{K-,p+) - '-bT{K-,p+) 



A^iB-^K-oj) = ^fBfKu{Xub2{K-,oj) + {Xu + X,) [63 (i^" , w) + br{K-,uj) 

A^iB-^TT^K*-) 
A^iB-^K-p^) 

A'^iB'^n^m*') 



2 

G_F 

2 

Gp 



fBUfK'\Xub2{K*',n'>) + (A„ + A,) 63(i^*~,7r°) 



62(i^-,P°) + (A„ + Ac)[63(i^-,P°) + 6r(^",P°) 



} 



/B/>/K.<!(A„ + Ae) 63(i^*°,r,W) 



■\bT{K*\r,^'^) + 'j^{bs{v^'\K*') - ^bTiv^'KK*") 



,\ ^*0, 1, 



A%B-^rj^')K*-) = ^/B/;,,/K.|A„[62(i^*-,r?W) + ^^b^iv^'^K*- 



(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 



(B15) 



+(A„ + Ae) [b:,{K*-,r,^'y) + bT{K*- ,r,^'^) + {b:,{r,^'\K*-) + bT{r,^'\K*-))\ | 



(3) Pure penguin processes: 

^''(S-^TT-T?*") = ^/B/./K*|A„62(i?*°,7r-) + (A„ + Ac)[63(:^*°,7r-) 

^«(S-^:^p-) = ^fBfKfp^Xub2{Tf,p-) + {Xu + Ae) [bs{T^,p-) + bT{Tf,P^ 

A-{B-^K-K*°) = ^fsfKfK^ [x'ub2{K*\ K-) + (A^ + A^) [b^{K*°, R-) + 6^ (if K-) 
A'^iB-^K^'K*-) = ^fsfKfK^ I^X'MK'',K*-) + (A^ + A'J [bs{K'',K*-) + bT{K\K*-) 

A''{B-^TT-(I)) = yl"(B°->7rV) = 0. 

^'^(B-^i^-,^) = ^fBfKf4>^Xub2{(t>,K-) + {Xu + Xc)[bz{(^,K-) + bT'{ct>,K-) 



] 



A'^(b'^K%) = ^fBf. 



>s:/<^|(A« 



+ Ac) bs{cl>,K'')-^br'{cl>,K") 



(B16) 
(B17) 
(B18) 

(B19) 

(B20) 
(B21) 

(B22) 
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